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Abstract 

Bovine heart eytoehrome c has been immobilized 
into the mesoporous siliea host material SBA-15 
in both its native folded and urea-unfolded state. 
The comparison of the two folding states’ behav¬ 
ior casts doubt on the commonly used explana¬ 
tion of cytochrome c adsorption, i. e. the electro¬ 
static interaction model. A detailed investigation 
of the protein binding as a function of pH and 
ionic strength of the buffer solution reveals the 
complex nature of the protein-silica interaction. 
Electrostatic interaction, van der Waais forces and 
entropic contributions by counterion release each 
contribute to adsorption on the silica pore walls. 

Introduction 

Over the course of the last two decades, protein 
adsorption in mesoporous materials has been in¬ 
tensively studied and plenty of technical applica¬ 
tions have been envisaged. Some proteins show 
an enhanced stability against denaturating condi¬ 
tions - chemical as well as thermal - and retain or 
even increase their electrochemical activity when 
immobilized in silica mesopores.^ More trivially, 
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since microorganisms like bacteria or fungi are far 
too large to penetrate mesoporous structures, en¬ 
capsulated proteins are well protected from bio¬ 
logical decomposition. This opens a wide field 
of biochemical applications that employ the enzy¬ 
matic activity of proteins under conditions which 
would otherwise destroy the enzymes. Other 
interesting applications arise from the fact that 
not all polypeptides adsorb equally well on all 
surfaccs.l^MI Thus, fractionation of complex pro¬ 
tein solutions and applications in chromatography 
should be feasible with customized mesoporous 
host materials. A third and rather promising field 
of interest is the use of porous materials as novel 
devices for controlled in-vivo drug release. Sev¬ 
eral researchers suggest utilizing this by loading 
a porous structure with drugs or enzymes and 
then injecting the loaded particles into living cells 
where the physiological conditions cause a release 
of the drugs into the cytosol. The feasibiliW of 
this application is demonstrated in literaturel^ and 
especially by Slowing et al.V^ The latter group 
loaded MCM-41 silica particles with a fluorescent 
protein which was subsequently released into the 
cytoplasm of human cervical cancer cells (HeLa- 
cells). 

As a general rule of thumb, protein adsorption 
appears to be non-specific and can be reversible 
under specific conditions. Interestingly enough. 
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despite the vast amount of researeh in this field, 
there is to date no eonsensus in literature what 
type of interaetion dominates the adsorption of 
biomoleeules on inorganie surfaees. It is the main 
motivation of this study to eompare two of the 
most eommon models of protein adsorption and to 
highlight their advantages and shorteomings. 

A widely used model system is the adsorption 
of eytoehrome c into the mesoporous siliea mate¬ 
rial SBA-15. There exists a vast amount of stud¬ 
ies eoneeming the adsorption of eytoehrome c on 
SBA-15, ineluding works by Deere)^, Hudsonl^, 
Miyahara^, Vinu^ and Zhang.^ Reviews were 
published by Hartmann^ and Zhao.® The behav¬ 
ior of this system is eommonly attributed to sim¬ 
ple eleetrostatie attraetion. Nevertheless, all these 
studies to date foeus on folded, i. e. native ey¬ 
toehrome c exelusively. Here, we propose a eom- 
prehensive analysis of both folded and unfolded 
eytoehrome c to gain new insights into the interae- 
tions whieh drive protein adsorption. 

Electrostatic Interaction Model 

The highest pore loadings, i. e. the amount of 
bound protein per gram of the siliea material, are 
often found elose to the isoeleetrie point where the 
overall eharge of the protein is zero. For exam¬ 
ple, Vinu et aZ.® studied the pH-dependent ad¬ 
sorption of horse heart eytoehrome c on SBA-15 
in 25 mM buffer solutions. They examined the pH- 
dependeney between pH 3 and pH 10.6. The high¬ 
est pore loading was observed at pH 9.6 whieh is 
only slightly below the pi of eytoehrome c. 

This is often interpreted in terms of a bal- 
aneing between an attraetive protein-wall interae¬ 
tion and protein-protein repulsion.®® The loss 
of eleetrostatie repulsion between the moleeules 
at their pi faeilitates the observed dense paek- 
ing of the adsorbing moleeules. While the over¬ 
all eharge of the protein vanishes at the pi, its 
surfaee still eontains patehes of positively and 
negatively eharged amino aeid residues. These 
eharged patehes drive the attraetion to the nega¬ 
tively eharged surfaee.®® At pH values far from 
the isoeleetrie point, the proteins will repell eaeh 
other and thus eause a less eompaet paeking den¬ 
sity on the adsorbing surfaee. 

The validity of this model was examined experi¬ 


mentally® and theoretieally®by studying the ad¬ 
sorption of lysozyme and a-laetalbumin on dif¬ 
ferently eharged surfaees. Taking into aeeount 
the ehemieal properties of the ionizable groups 
and the orientation of the adsorbed protein with 
respeet to the surfaee, the results of the exper¬ 
imental study were reprodueed in a quantitative 
manner by just eonsidering eleetrostatie interae¬ 
tion. Any other effeets like dispersion forees, hy- 
drophobie interaetions, eonformational ehanges or 
any other type of interaetion were not required to 
reproduee the experimental findings. One eould 
thus assume that eleetrostatie interaetion between 
eharged patehes, net eharge and surfaee eharges is 
the dominant meehanism behind protein adsorp¬ 
tion. Nevertheless, this interpretation is laeking 
at one point. The experiments employ surfaees at 
least partially eovered with a tethered PEG spaeer 
layer with a thiekness exeeeding the Debye length 
of the buffer. Contrary to the experimental reality, 
the theoretieal deseription eonsideres an uneov- 
ered, blank surfaee. Thus the exaet nature of the 
protein-surfaee interaetion still remains unelear. 

However, eleetrostatie interaetion still seems to 
be exeeptionally suitable to explain the behavior 
of folded eytoehrome c on negatively eharged sil¬ 
iea.® 

Counterion Release Model 

When a eharged plate is immersed into a solution 
eontaining positive and negative ions the eoneen- 
trations of these ions will ehange in the vieinity 
of the surfaee. Ions with the same sign as the 
surfaee eharges (eoions) will be repelled into the 
bulk while those with a different sign (eounteri- 
ons) will be drawn towards the surfaee. This leads 
to the formation of a shielding layer. For low ionie 
strength eleetrolytes and highly eharged surfaees, 
the eounterion eoneentration in the shielding layer 
strongly surpasses the bulk eoneentration. This 
leads to a eonsiderable entropy loss for the bound 
ions. When a seeond plate with opposite eharge is 
brought in eontaet with the first plate, their eharges 
mutually shield eaeh other and the eounterions are 
released into the bulk, sinee the shielding layer is 
no longer needed to ensure eleetroneutrality. This 
is aeeompanied by an entropy gain whieh eauses 
an attraetive foree between the plates that adds to 
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the mere eoulombie attraetion. Explieit ealeula- 
tions^^of the resulting forees highlight the impor- 
tanee of this entropy-driven interaetion espeeially 
for low ionie strength solutions. 

This model ean be easily extended to protein ad¬ 
sorption, at least in a qualitative manner: Both the 
siliea surfaee and the eharged patehes on the pro¬ 
tein are eovered with a shielding ion layer whieh 
desolves upon adsorption. For low ionie strength 
eleetrolytes, the entropie attraetion of the eoun- 
terion release should thus play an important part. 
Furthermore, if adsorption is dominated by eoun- 
terion release, we ean expeet the binding affinity to 
drop signifieantly if the ionie strength inereases. 
This behavior has been observed experimentally 
for a multitude of different protein-surfaee eom- 
binations like e.g. hemoglobin on ola)^, ey- 
toehrome c on eyano-funetionalized SBA-IS^^, ey- 
toehrome c on fused siliea^ and lysozyme in 
eharged mierogels.l^ 

Theoretieal work eonfirms that this eounte- 
rion release meehanism ean indeed explain the 
experimentally observed Fangmuir-type binding 
isotherms .1^ 

Experimental 

Cytochrome c 

Bovine heart eytoehrome c was purehased from 
Sigma Aldrieh, eatalog number C2037, and used 
as-reeeived without further treatment or purifiea- 
tion. It eonsists of a single polypeptide ehain 
of 104 amino aeid residues whieh are eovalently 
bound via two eystein residues to a eentral heme 
eomplex. Its geometrieal dimensions in the na¬ 
tive, folded state have been reported^ as 26 x 
“ 3 

32 X 33 A . Immersion into 8 M urea solutions dis¬ 
rupts the hydrogen bonds whieh stabilize the ey¬ 
toehrome c’s tertiary strueture and leads to unfold¬ 
ing of the moleeule.^ Small angle x-ray data^^ 
of the unfolded eytoehrome c at pH 7 reveal a 
struetural transition from the almost spherieal eon- 
formation of the folded protein (semi-major axis 
18 A, semi-minor axis 18 A, radius of gyration 
Rg= 12.8 A) to an eeeentrie ellipsoid shape (semi¬ 
major axis 65 A, semi-minor axis 9 A at 8 M urea 
and Rg = 29.7 A at 10 M urea). 


Preparation of Mesoporous SBA-15 Sil¬ 
ica 

The synthesis of hexagonally ordered mesoporous 
SBA-15 was first reported by Zhao et al.^ The 
samples used in this thesis were prepared aeeord- 
ing to the following proeedure: We mix 4 g of the 
tri-bloek eo-polymer PEO 20 -PPO 70 -PEO 20 with 
129.6 g water and 19.3 ml HCl (37 %). Due to its 
amphiphilie nature, the polymer forms an ordered 
phase of mieellar struetures when mixed with wa¬ 
ter. Vigorous stirring at 350 rpm for four hours 
is needed to ensure a homogenous emulsion. The 
mixture is kept in an oil bath at 55 °C during this 
proeess. We then add 8.65 g tetraethylorthosili- 
eate (TEOS) and stir the system for another 20 
hours. We subsequently inerease the temperature 
to 85 °C and let the mixture rest for another 22 
hours without stirring. During this time, the silieon 
from the TEOS leads to an aeeumulation of siliea 
around the polymer mieelles. These aggregates 
preeipitate as a fine-grained powder. Caleination 
of the repeatedly rinsed powder at 500 °C finally 
removes the polymer while preserving a negative 
of the mieellar strueture in the siliea grains. The 
porous siliea powder ean now be used without fur¬ 
ther treatment or purifioation. Small angle x-ray 
diffraetion patterns were reeorded at DESY, Ham¬ 
burg. Five Bragg peaks were observed, eonfirming 
a hexagonal arrangement of linear mesopore s with 
a lattiee parameter of a/, = 10.71 ± 0.08 

Nitrogen Sorption Isotherms 

Nitrogen sorption isotherms were eondueted by 
eontrolled filling (adsorption) and evaeuation (des¬ 
orption) of the sample via a eustom-made gas han¬ 
dling system. The main part of the gas handling 
system was kept at room temperature while the 
sample eell was eooled to a well-known referenee 
temperature. This was aehieved by either using 
a elosed-eyele helium refrigerator (Eeybold RGD 
510 Cryostat with RW 2 Compressor Unit) or by 
simply immersing the sample eell into a dewar 
vessel filled with liquid nitrogen. The pressure 
relaxation inside the system was measured using 
a Baratron Capaeitanee Manometer (MKS Intru- 
ments) with 1000 torr full seale. 
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Cytochrome c Adsorption to SBA-15 

Eppendorf Safe-Lock tubes were filled with small 
amounts of SBA-15. The mass of these SBA- 
15 samples (usually 4.0 — 9.0mg) was mea¬ 
sured with an analytical balance (Sartorius type 
1801). To provide for proper stirring during 
the experiment, we equipped each tube with 
a small stirring bar. Prior to the insertion 
the bars were rinsed with acetone and subse¬ 
quently air-dried at 60°C. Cytochrome c solutions 
were prepared with the following concentrations: 
2000 mg/1, 1500 mg/1, 1250 mg/1, 1000 mg/1, 
750 mg/1, 500 mg/1, 375 mg/1, 250 mg/1, 125 mg/1 
and 42 mg/1. For each concentration two tubes 
were filled with 200 /il of protein solution per 
1 mg of SBA-15. To ensure a fine dispersion of 
the silica grains, the samples were immersed into 
an ultrasonic bath for 10 minutes. To remove any 
impurities the solutions were filtered with dispos¬ 
able 0.8 jum syringe filters prior to the adsorption 
experiments. 

During the adsorption process, the samples were 
kept in a heated water bath at SEC and stirred 
at 350 rpm using a IKA RCT basic safety control 
magnetic stirrer. According to similar procedures 
reported in literature,^^ we let the system equi¬ 
librate for 5 days. Estimations of the adsorption 
kinetics indicate that this time is sufficient to reach 
equilibrium. 

The amount of cytochrome c bound to the sil¬ 
ica was calculated from the difference in the initial 
and final protein concentration. 



Figure 1: Top: Absorbance spectra of ferric cy¬ 
tochrome c in pure acidic buffer (red line) and in 8 M 
urea acidic buffer (black line). Bottom: Spectra of fer¬ 
rous (black) and ferric (red) cytochrome c. Ferrous cy¬ 
tochrome c shows two distinct peaks at 520 and 550 nm. 

Ionic Strength dependent Adsorption 

Protein solutions with a fixed concentration of 1 g/1 
were prepared using buffers with an initial ionic 
strength of 10 mM. We adjusted the ionic strength 
of each sample by mixing approriate amounts 
of pure and salinated buffers containing 1.11 M 
NaCl. The addition of urea and salt altered the 
initial pH values of the buffers. However, these 
changes were found to be marginal and are there¬ 
fore considered negligible. There is, however, 
one exception: The acidic buffer experiences quite 
substantial alterations of its pH value when sali¬ 
nated in the absence of urea. The implications 
of this changes will be discussed where necessary. 
Apart from the different preparations of the solu¬ 
tions, the ionic strength dependent measurements 
were carried out according to the same procedure 
as the adsorption isotherms mentioned above. 
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Concentration Measurement via UV- 
Vis Spectroscopy 

The supemate concentration after the adsorption 
was determined photometrically with a Hitachi U- 
3501 spectrophotometer and disposable UV-Vis 
cuvettes {Plastibrand PMMA 2.5 ml macro cu¬ 
vettes). Absorbance spectra of cytochrome c, 
i. e. the molar absorbtivity as a function of the 
wavelength a = a(A), are shown in figure 
The spectra are altered by chemical and structural 
transitions of the protein.l^ We used the strong 
Soret band absorption arround 410 nm to simul¬ 
taneously check the concentration and the confor¬ 
mational state of the protein. Additionally, we 
used the Q band region between 500 and 600 nm 
to monitor the oxidation state of the protein: Sam¬ 
ples containing reduced cytochrome c show two 
distinct peaks in this region. Since proper un¬ 
folding can only be achieved with oxidized cy¬ 
tochrome cm, any samples showing signs of these 
peaks were rejected from further evaluation. In ac¬ 
cordance with literaturel^l^, complete unfolding 
of the oxidized protein was assumed in urea solu¬ 
tions without further experimental confirmation. 

Buffer Solutions 

All buffers were prepared in-house according to 
the following recipes. PH-dependent measure¬ 
ments used buffers with pH 3.0 (0.018 M cit¬ 
ric acid and 0.033 M trisodium citrate), pH 6.0 
(0.025 M monopotassium phospate and 0.002 M 
sodium hydroxide), pH 9.0 (0.018 M glycin, 
0.017 M sodium chloride and 0.03 M sodium hy¬ 
droxide). The high urea content needed to unfold 
the proteins considerably altered the pH value 
of the buffers to 4.42, 6.40 and 9.69, respec¬ 
tively. To ensure comparability between the mea¬ 
surements of folded and unfolded cytochrome c, 
new buffers were prepared for the measurements 
without added urea as follows: pH 4.4 (0.01 M 
trisodium citrate and 0.011 M citric acid) and pH 
9.7 (0.0135 M glycin, 0.014 M sodium chloride 
and 0.073 M sodium hydroxide). The pH change 
in the near-neutral buffer was found small enough 
to use the same buffer for the measurements with 
and without urea since both the initial and the 
altered pH were comfortably in the region were 


protein and surface bear opposite charge signs. 
Ionic strength dependent measurements used the 
following buffers: pH 3.0 (116.37 g 0.01 M cit¬ 
ric acid and 4.0 ml 0.01 M trisodium citrate), 
pH 3.8 (200.02 g 0.01 M acetic acid and 15 ml 
0.01 M sodium acetate), pH 4.5 (100.23 g 0.01 M 
trisodium citrate and 120 ml 0.01 M citric acid), 
pH 6 (100.14 g 0.01 M monopotassium phospate 
and 32.5 ml 0.01 M sodium hydroxide), pH 7.3 
(50.40 g 0.01 M sodium hydroxide and 41.7 ml 
0.05 M monopotassium phospate), pH 8.5 (83 ml 
0.01 M sodium tetraborate and 13.1ml 0.01 M 
hydrochloric acid) and pH 10.6 (100 g 0.01 M 
sodium bicarbonate and 71 ml 0.01 M sodium hy¬ 
droxide). 
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Results 

SBA-15 Characterization 



Figure 2: Scanning electron micrograph of SBA-15 
powder. The contrast of all pictures was enhanced with 
standard image processing software (GNU Image Ma¬ 
nipulation Program GIMP 2.6.5). Top: Aggregate of 
individual grains. Middle: Lateral view of a single 
grain. Bottom: The front edge of a silica grain. 

Figure shows scanning electron micrographs 
of the SBA-15 powder used in this study. The 
powder consists of needle-like aggregates of 
micrometer-sized silica grains. A zoom in on 
the edge of such a single grain reveals the porous 
structure. The picture in the bottom panel shows a 
micrograph of the front edge of a grain. The pore 
openings are visible as dark spots on the brighter 


silica surface. The middle panel presents a lateral 
view of the grain surface where the pores can be 
identified as dark grooves. Their high aspect ratio 
is evident in this micrograph. 

Nitrogen sorption measurements were per¬ 
formed as outlined above. A sorption isotherm 
is the plot of the total amount of adsorbed nitrogen 
n = Y^Ansamp versus the corresponding relaxation 
pressure prelax- For the sake of simplicity and 
comparability, we usually use the dimensionless 
quantities filling fraction / = h/hq and reduced 
vapor pressure P = p/pQ. Here, hq is the amount 
of nitrogen needed for a complete filling of the 
sample’s micro- and mesoporous structures. As 
can be seen in the sorption isotherm shown in fig¬ 
ure there is some condensation beyond hq. This 
can be attributed to highly irregular macroporous 
structures, po denotes the bulk vapor pressure at 
the sample temperature. 



Figure 3: Experimental nitrogen adsorption (trian¬ 
gles) and desorption (circles) isotherm and theoretical 
isotherms according to the Saam and Cole Theory for 
monodisperse pores (dash-dotted line), Gaussian pore 
size distribution (dashed line) and the bimodal pore size 
distribution (solid line) shown in the inset. This figure 
is taken from a previous study 

At very low pressures the interaction between 
the gas molecules and the naked surface causes 
the gas to condensate. The large initial slope of 
the isotherm reflects this phenomenon. The kink 
in the isotherm at a reduced pressure of approx¬ 
imately 0.05 marks the point were a complete 
monolayer has formed. Further increasing of the 
pressure leads to a growth of the condensed layer’s 
thickness. This region which spans up to reduced 
pressures of 0.6 is called the reversible multilayer 
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regime. Upon reaehing a eritieal thiekness the 
liquid film beeomes unstable. Further adsorption 
of gas moleeules then leads to the spontaneous 
formation of capillary bridges, i. e. small liquid 
droplets whieh span aeross the entire pore diam¬ 
eter. In figure lathis transition ean be seen for p/po 
between 0.6 and 0.7 as a strong inerease in the 
isotherm’s slope. This eondensation is a first or¬ 
der phase transition. Adding more gas moleeules 
to the system will eause the droplets to grow longi¬ 
tudinal to the pore axis while the pressure remains 
unaltered. Note that there is a strong hysteresis 
in the pressure assoeiated with eapillary eonden¬ 
sation between adsorption and desorption. 

We analyze the reeorded isotherms using the 
mean-field model by Saam and Cole.^l^ Con¬ 
trary to mueh simpler models like the BET or the 
BJH approaeh, this theory explieitely ineludes the 
van der Waals interaetion between eondensate and 
pore wall. 

The green dash-dotted line in figure eorre- 
sponds to the hypothetieal isotherm of a sample 
with a monodisperse mesopore radius of rp = 
3.3nm, while the blaek dashed eurve ineludes a 
Gaussian distribution of pore sizes. We ealeulated 
150 individual Saam and Cole isotherms, eaeh 
with a different pore radius but otherwise identieal 
parameters. The final isotherm is a sum of these 
isotherms, weighted aeeording to a single Gaus¬ 
sian distribution and normalized to a filling frae- 
tion f = I dX p/pq = 1. It reproduees the exper¬ 
imental hysteresis loop rather well, but severely 
underestimates the multilayer region. This ean 
be fixed by using a bimodal pore size distribu¬ 
tion using two Gaussian distributions. Refering 
to the SAXS data by Imperor-Clere et al^, we 
attribute those two peaks to the hexagonal meso¬ 
pore array and a mieroporous layer surrounding 
the individual pores. We also added a lower eut- 
off of 0.5 nm to the pore sizes allowed in the eal- 
eulation. This eut-off was ehosen for numeri- 
eal reasons, but it is also reasonable to assume 
a physieal eut-off for the SBA-15’s mieroporos- 
ity: Sinee the mieropores presumably stem from 
individual hydrophilie polymer segments, no mi- 
eropore should be smaller than the baekbone of 
the polymer ehains. The samples used in this 
study show a broad peak at r„,/cro = 0.75 nm with 

micro = 78% and a narrower mesopore peak at 
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3.3 nm with Omeso = 6.5%. 


PH-Dependent Cytochrome c Adsorp¬ 
tion 


The results of the protein adsorption experiments 
are shown in figure The experiments were 
done at three different pH values eorresponding 
to three fundamentally disparate eleetroehemieal 
eonditions. The red line and symbols eorrespond 
to a pH of 9.7, elose to the isoeleetrie point of the 
eytoehrome c. The protein is virtually uneharged 
under these eonditions, while the siliea exhibits 
strong negative surfaee eharges. The green line 
and symbols were measured under near-neutral 
eonditions where both the protein and the surfaee 
are eharged with opposite signs, resulting in mu¬ 
tual attraetion. Finally, the blue lines and sym¬ 
bols represent the measurements near the SBA- 
15’s isoeleetrie point. The protein has a strong 
positive overall eharge at this pH, while the surfaee 
is mostly neutral (The isoeleetrie point of SBA-15 
has been reported as pH3.8.l^. The solid lines in 
the figure are Langmuir-type fits defined by 
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This equation relates the amount of adsorbed 
protein n^ds to the residual supernate eoneentration 
Cs, using the amount of protein ni needed to form 
a eomplete monolayer eovering the sample surfaee 
and an interaetion parameter a as fit parameters. 
This behavior is well-known: Isotherms similar to 
the ones presented here were already presented by 
Vinu et al.^ 
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Figure 4: Cytochrome c adsorption isotherms for pure 
buffer (top) and for 8M urea buffer (bottom). Solid 
lines are Langmuir-type fits to the experimental data. 
The fit parameters are listed in table [T] The isotherms 
were recorded pH 9.7 (red squares), pH 6.4 (green cir¬ 
cles) and pH 4.4 (blue triangles). 

For both folding states of cytochrome c the pore 
loadings, i. e. the amount of adsorbed protein, de¬ 
crease with decreasing pH value. This decrease 
is quite moderate for unfolded cytochrome c but 
severe for the native, folded type. While the na¬ 
tive cytochrome c adsorbs up to 25.9±4.3 /imol 
per gram of SBA-15 near its pi and still up to 
22.9±7.4/imoFg at pH 6.4, the adsorption is al¬ 
most negligible on the weakly charged surface 
at pH4.4 with a mere 3.4T0.4 /imol/g saturation 
value and an interaction parameter a = 0.023 ± 
0.006 l//imol, almost 20 times smaller as for the 
alkaline buffer. Unfolded cytochrome c, on the 
other hand, seem less susceptible to changes in the 
buffer pH. Both the saturation value of the pore 
loading and the interaction parameter at pH 4.4 are 
still about half of the corresponding values at the 
protein’s pi. 


To evaluate what percentage of the available 
pore space is actually occupied by adsorbed 
molecules, the volume and surface of the meso- 
pores were extracted from the nitrogen sorp¬ 
tion data and the calculated pore size distribu¬ 
tion. The samples have a mesopore volume of 
ymeso = 0.520cm^/g and a specific surface area 

Ameso = 317m^/g. The geometrical dimensions!^ 

° 3 

of folded cytochrome c (26 x 32 x 33A ) yield a 
volume of Vy = 14.4nm^ per molecule. Assuming 
a hexagonal packing of the individual molecules 
in a surface-covering monolayer, each protein 
takes up 11.2nm^ of the surface area. A complete 
and perfect monolayer would thus correspond to 
1.48 ■ 10^^^ mol/cm^, which is in accordance with 
the literature.!^ With these data, the highest pore 
loading of our experiments (25.9 /imol/g) implies 
a packing density of 8.1 -lO^^^mol/cm^ taking up 
43% of the mesopore volume. This is only about 
half of the expected value, but since the highly 
curved and probably corrugated pore surface is 
unlikely to allow a perfect hexagonal protein as¬ 
sembly this deviation is within reasonable lim¬ 
its. We also expect the small size of the pores 
to impose a steric hindrance to perfect mono- 
layer adsorption. It is further possible that the 
protein experiences conformational changes upon 
adsorption which increase the area that is occu¬ 
pied by each single molecule. Our value for the 
percentage of the pore volume which is filled upon 
adsorption coincides well with previous findings 
by Miyahara et al.^ They found that adsorption 
of cytochrome c takes up approximately half the 
pore volume of SBA-15. Given the respective 
radii of the protein and the pores, they calcu¬ 
late that a close packing of spherical molecules 
inside a cylinder takes up 48% of the available 
cylinder volume. We therefore assume that in the 
case of adsorption of the native protein at pH 9.7 
the pores of our sample are completely stuffed 
with cytochrome c. At pH 6.4 the surface and 
volume packing densities are 7.1 -lO^^^mol/cm^ 
and 38%, respectively. At pH 4.4, we observe a 
surface coverage of 1.06-lO^^^mol/cm^ and a 
volume packing density of 5.6%. 

We obtain the respective values for the unfolded 
protein using the x-ray scattering data of Hsu et 
al^ In 8M urea solutions, each molecule has 
a volume of Vu = 22nm^ and takes up 23.4nm^ 
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when adsorbed on a surfaee. The theoretieal 
monolayer would thus eontain 7 ■ lO^^^mol/em^. 
Again, the experimental data are about half of 
this value. At pH 9.7 we observe a eover- 
age 4.0-lO^^^mol/em^ and a paeking density of 
33%. At pH6.4 we observe 2.8-lO^^^mol/em^ 
and 23% paeking density while pH 4.4 yields 
2.5 -lO^^^mol/em^ and 20%, respeetively. 

Influence of the Ionic Strength 

Figure depiets the adsorbed amount of ey- 
toehrome c as a funetion of the ionie strength. 
The data refer to the aetual number of proteins 
adsorbed from an 1 g/1 solution and must not be 
eonfused with the monolayer eoverage ni from 
the Langmuir-type fits. However, 1 g/1 equals an 
initial eoneentration of approximately 80 /imol/1. 
In the limit of low adsorption where the initial 
and residual supernate eoneentration differ only 
slightly this should be suffieient to ensure that 
the Langmuir isotherm is almost saturated. Apart 
from the folded eytoehrome c samples in alkaline 
buffer, it is therefore not too far-fetehed to use 
the measured adsorbed amount as an estimate of 
ni. The pH values differ slightly from the ones 
used in the previous experiments but still eorre- 
spond to the two distinet isoeleetrie points and 
the near-neutral eonditions. Blaek symbols show 
the measurements with pH 10 buffer, red symbols 
used pH 6.5 buffer and the blue symbols eorre- 
spond to pH 3.8 buffer. The folded protein data 
for the alkaline buffer exhibit a maximum at an 
ionie strength 0.2 M/1 and a deerease to approxi¬ 
mately 10/imol/g at 1 M/1 NaCl. The data from 
the pH 6.5 buffer do not show this maximum but 
also a monotonie deerease to 7 /imol/g at 0.4 M/1 
and beeome only weakly suseeptible to further 
ehanges in the salinity. The measurements us¬ 
ing the aeidie buffer showed a quite different be¬ 
havior. Starting at roughly 2 jumol/g at low salt 
eontent, the adsorbed amout of protein inereases 
almost linearly to 5 /rmol/g at 1 M/1. The ad- 
sorbtion of unfolded eytoehrome c is mueh more 
strongly influeneed by salinity. Adsorption from 
alkaline buffers deereases strongly until saturating 
at about 2.5 /imoFg above 0.7 M/1 NaCl. Adsorp¬ 
tion from near-neutral buffers even beeomes al¬ 
most negligible for 0.3 M/1 and higher. Again, the 


aeidie buffers show a different behavior. Both for 
very high and low salinities, the adsorption ranges 
around 2.0-2.5 /rmol/g, but undergoes a minimum 
of approximately 1 /rmol/g at 0.4 M/1. 




ionic strength [M/I] 


Figure 5: Equilibrium value of cytochrome c adsorp¬ 
tion from a 1 g/1 solution as a function of ionic strength 
for pure buffer (top) and for 8 M urea buffer (bottom) at 
pH 10.0 (black squares), pH 6.5 (red circles) and pH 3.8 
(blue triangles). 


Discussion 

Electrostatic Interpretation 

As expected from previous studies,!^^ the ad¬ 
sorption of folded cytochrome c can be readily 
explained with the electrostatic interaction model. 
The highest pore loading is found at the isoelectric 
point of cytochrome c, where the protein-protein 
repulsion vanishes and the protein-wall attraction 
facilitates high packing densities. At intermedi- 
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ate pH, the packing density decreases due to the 
rise of electrostatic protein-protein repulsion. Fi¬ 
nally, at the silica’s pi, protein adsorption is largely 
supressed due to vanishing silica surface charges. 
This is not only reflected in the very small satura¬ 
tion value of the Langmuir isotherm nQ pH= 4 A = 
3.37 /rmol/g, but also in the severe decrease of the 
interaction parameter a, which drops by almost a 
factor of 20 from 0.40±0.191//imol at pH 9.7 to 
0.023 ± 0.0061/Atmol at pH 4.4. 

This picture also applies to the adsorption of the 
unfolded protein. Again, the highest pore load¬ 
ing is observed at the protein’s pi, while the bind¬ 
ing affinity is rather low at pH 4.4. The maxi¬ 
mum pore loadings of unfolded protein are ap¬ 
proximately half the value of the folded type. This 
fits well with the surface area occupied by a sin¬ 
gle molecule, which increases from 11.2nm^ to 
23.4 nm^ upon unfolding. 

There is, however, one puzzling point: The 
isotherm of unfolded cytochrome c in the acidic 
buffer does not experience the severe drop in both 
its parameters compared to the folded protein. 
Electrostatic interaction alone fails to explain this 
behavior. We found it plausible to assume that van 
der Waals forces cause this different behavior. At 
the silica’s pi they probably dominate the protein 
binding due to the absence of electrostatic inter¬ 
action. The folded protein has a relatively rigid 
shape. While some data suggests small conforma¬ 
tional changes upon adsorptionl^, these changes 
are rather minor. As a consequence, the folded 
protein has a rather small contact area with the sil¬ 
ica and can not exploit these interactions very ef¬ 
fectively. The unfolded protein, on the other hand, 
has a much higher contact area with the silica and 
is probably rather flexible.l^It should therefore be 
much more capable to exploit the van der Waals in¬ 
teraction and still be able to adsorb to a relatively 
high extend even in the absence of electrostatic at¬ 
traction. 

To test this hypothesis, we examined the influ¬ 
ence of the buffer’s salinity on the protein adsorp¬ 
tion. If it is true that the unfolded protein’s adsorp¬ 
tion behavior is influenced by van der Waals forces 
and the folded protein’s behavior is solely due 
to electrostatic interaction, changes in the salinity 
should affect the unfolded samples less than the 
folded ones. A single glance at figure [^reveals that 


the opposite is true: The adsorption of unfolded 
cytochrome c is highly susceptible to changes in 
the ionic strength and becomes almost negligible 
above 0.7 M/1 NaCl. Folded cytochrome c, on the 
other hand, adsorbs still quite well even at high 
salinity, especially at the protein’s pi. 

Interfacial charge regulation has to be taken into 
account when discussing electrostatically driven 
adsorption. As discribed in detail by Hartvig et 
al^, the enrichment of counterions over a charged 
surface leads to a alteration of the pH-value in the 
Debye layer. This altered pH-value will affect the 
local charges on the protein and might lead to a 
behavior which differs from what would be ex¬ 
pected from the bulk solution’s properties. How¬ 
ever, a quantitative inclusion of this mechanism 
into our interpretation encounters two problems. 
First, the calculation of the surface pH requires 
knowledge of the surface potential, which is in¬ 
fluenced not only by the protonation of the sil¬ 
ica’s silanol groups, but also on the amount and 
charge of the adsorbed proteins. Second, even 
for a known alteration of the surface pH the im¬ 
pact on the protein adsorption can only be calcu¬ 
lated with precise knowledge of the pH-dependent 
charge profile of the protein, its structure, the spa¬ 
tial distribution of charged surface patches and the 
orientation of the molecule upon adsorption. To 
the best of our knowledge, these data do not ex¬ 
ist for unfolded cytochrome c. We can, however, 
present some qualitative considerations. At the 
acidic pH, the charge of the silica surface vanishes 
and thus charge regulation should not play a ma¬ 
jor role. At higher pH, the adsorption of the (pre¬ 
sumably positively charged) proteins will screen 
the negative surface charge to a certain extend and 
suppress charge regulation. The fact that the Lang¬ 
muir isotherms for cytochrome c presented here 
and in literature can be readily explained with¬ 
out assuming charge regulation suggest the it plays 
only a negligible role in the model system used 
here. 

Counterion Release Interpretation 

The surface of a protein is composed of posi¬ 
tively and negatively charged patches of different 
size as well as polar and hydrophobic regions.!^ 
The counterion release mechanism can be used 
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to explain protein adsorption to eharged surfaees 
simply by eonsidering the loeal shieldiim of the 
eharged surfaee patehes. Beeker et al^ argue 
that eounterion release is an effeetive meehanism 
as long as the ionie strength of the bulk solu¬ 
tion is smaller than the surfaee eharge eoneentra- 
tion of the eharged patehes. As an example they 
mention /3-laotoglobulin whieh has 5 eharges in 
a pateh of 10 nm^, eorresponding to an ion eon- 
eentration of 0.3 M, assuming that the shielding 
ions are eonfined within one Debye length from 
the surfaee. Thus, any attraetion from eounterion 
release should vanish for higher ionie strength. Of 
eourse the exaet values might be different for ey- 
toehrome c, but without preeise knowledge of the 
eytoehrome’s surfaee eharges at the various ehem- 
ieal eonditions examined we might still use the 
given example as a reasonable, semi-quantitative 
referenee. 

With this value in mind, we ean interpret the 
data from figure If we negleet the data for the 
aeidie buffers we notiee the general trend of de- 
ereasing adsorption for inereasing salinity. This 
deerease seems not very strong for folded protein 
in the alkaline buffer. But note that the starting 
eoneentrations were too low to saturate the ad¬ 
sorption, so the data might underestimate the true 
extent of this deerease. Adsorption of folded ey- 
toehrome c at pH 6.5 strongly deereases for in- 
ereasing ionie strength and seems to saturate above 
0.4 M. This eoineides well with the estimated ion 
eoneentration limit of Beeker et al.. Any adsorp¬ 
tion at higher ionie strength is likely due to van 
der Waals forees. Unfolded eytoehrome c be¬ 
haves quite similarly: The eurve for pH 10 satu¬ 
rates somewhere between 0.5 M and 0.7 M while 
adsorption is merely over the deteetion limit for 
pH 6.5 above 0.2 M. These data indieate that eoun¬ 
terion release is indeed the driving meehanism be¬ 
hind eytoehrome c adsorption. 

However, there are still some features whieh 
the eounterion release model fails to aeeount for. 
The high pore loadings of folded eytoehrome c at 
pH 10 ean only be explained if we assume highly 
eharged, small surfaee patehes. The inerease of 
the pore loading for both eonformational states 
in aeidie buffers above 0.4 M also remains unex¬ 
plained. Coneerning the data from figure the 
absenee of the sudden drop in unfolded protein 


pore loading at pH 4.4 is still unaeeounted for and 
the higher pore loadings at higher pH values ean 
only be explained if we make the ad hoc assump¬ 
tion that the eharge density of the patehes inereases 
upon approaehing the isoeleetrie point. Not only 
would this assumption be overly tentative, it is 
also highly eounterintuitive given the faet that the 
overall eharge of the protein vanishes at the pi. 
We therefore eonelude that the eounterion release 
meehanism alone is also insuffieient to explain ey¬ 
toehrome c adsorption. 

Complete Interpretation 

None of the simple models presented above is suf- 
fieient to exelusively aeeount for the observed pro¬ 
tein adsorption behavior. The driving foree behind 
eytoehrome c immobili z ation on siliea is probably 
a rather eomplex eompilation of several interae- 
tions. We propose a eomprehensive explanation of 
the experimental data using eleetrostatie and eoun¬ 
terion release interaetion as well as van der Waals 
forees. 

The eleetrostatie interaetion ean aeeount for the 
high pore loadings elose to the proteins isoelee¬ 
trie point, whieh are rather puzzling if we eon- 
eem only the eounterion release. It furthermore 
explains the small maximum of the ionie strength 
dependent pore loading for folded eytoehrome c at 
pH 10. Even though the adsorption is unsatu¬ 
rated for these samples, the maximum at 0.2 M 
NaCl is unambiguous and indieates a higher bind¬ 
ing affinity. While residual interprotein repul¬ 
sion will probably play a role at very low ionie 
strength it is likely shielded at 0.2 M, faeilitating 
a denser paeking. Both the eounterion release and 
the eleetrostatie interaetion will suffer from in¬ 
ereasing ionie strength. The first will do so due 
to a deereasing entropy gain while the latter will 
simply be shielded by the solution’s high ion eon- 
eentration and thus shorter Debye length. This be¬ 
havior is very well refleeted in the ionie strength 
dependent measurements of figure The esti¬ 
mation of Beeker et al^ that the eounterion re¬ 
lease meehanism should be ineffeetive at salini¬ 
ties above 0.3 M eoineides well with our data, as 
outlined above. It is still not entirely elear why 
the unfolded protein is more suseeptible to salin¬ 
ity ehanges. Without preeise knowledge of the un- 
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folded protein’s strueture and orientation, we ean 
not rule out interfaeial eharge regulation as a pos¬ 
sible explanation. But given the eonsiderartions 
outlined above, we propose an interpretation us¬ 
ing the eounterion release meehanism: Sinee the 
unfolded moleeules have a mueh higher surfaee 
area it is not very far-fetehed to assume that their 
eharged amino aeid residues are distributed over a 
larger area. This would reduee the ion eoneentra- 
tion needed to shield the patehes and reduee the 
ionie strength at whieh the eounterion release be- 
eomes effeetless. This redueed effeetiveness of the 
eounterion release might also be refleeted in the 
Langmuir interaetion parameters a. With the ex- 
eeption of the native sample at pH 4.4, the values 
for a at a given pH are eonsiderably smaller for the 
unfolded samples. Sinee a represents the interae¬ 
tion strength, this is eonsistent with the assumption 
of a lowered eharge density in the unfolded state. 

At the siliea’s isoeleetrie point both eleetrostatie 
and eounterion release interaetion will eease to 
work. Any residual adsorption will therefore be 
dominated by van der Waals forees. Unfolded ey- 
toehrome c has a higher surfaee area and is proba¬ 
bly more flexible than its native eounterpart. This 
enables the protein to squeeze tightly to the surfaee 
and exploit the van der Waals interaetion rather ef- 
feetively. It ean therefore still adsorb quite well in 
the aeidie buffers whereas the folded protein ex¬ 
hibits a very low binding affinity. 

The peeuliar behavior of the ionie strength- 
dependent measurements for aeidie buffers is 
likely an artefaet of the buffers used. While the 
pure buffer with pH 3.86 is very elose to the sil¬ 
iea’s isoeleetrie point, adding 1 M NaCl severely 
reduees the pH to 3.2. Thus, inereasing salin¬ 
ity will lower the buffer pH below the siliea’s 
pi. The pore surfaee will exhibit positive eharges 
and allow for eleetrostatie binding with negative 
patehes on the protein, leading to the almost linear 
pore loading inerease oberserved for the blue data 
points in the upper panel of figure Additionally, 
higher salinities will enhanee the pore loading by 
shielding the interprotein repulsion. This behavior 
is slightly altered when urea is added. The pure 
buffer with urea has a pH of 3.96. At this pH, the 
siliea is still slightly eharged, allowing for some 
residual eleetrostatie binding. Adding 1 M NaCl 
alters the pH to 3.7, whieh is below the siliea’s 


pi. Again, a small amount of surfaee eharges will 
be present, although with different sign. At some 
salinity in between, the pi will be met exaetly, re¬ 
sulting in almost negligible binding. 

A last aspeet whieh remains unaeeounted for is 
the high pore loading of folded eytoehrome c at 
high salinity elose to its pi. The pore loading is 
still higher than 10 /imol/g, about twiee the value 
observed for the folded protein in the other buffers 
and about four times that of the unfolded samples. 
Interprotein repulsion should be shielded anyway 
at this ionie strength, so the differenee between the 
alkaline and near-neutral measurements ean not be 
explained within the seope of eleetrostatie inter¬ 
aetion. Van der Waals forees are not aeeountable 
for this, either, sinee they are independent of the 
ehemieal eonditions. And without preeise knowl¬ 
edge of the surfaee eharge distribution under these 
very speeifie ehemieal eonditions it is impossible 
to rule out eontributions from eounterion release 
or eharge regulation. Note that we assume that 
the protein strueture is identieal for pH 6.4 and 
pH 9.7. However, eytoehrome c undergoes the al¬ 
kaline transition at the latter pH. This transition is 
marked by the replaeement of the Met-80 ligand of 
the heme group by a lysine and results in a strue¬ 
ture whieh is still folded but more flexible than the 
native one.® Analogous to the explanation used 
for unfolded eytoehrome c this enhaneed flexibil¬ 
ity might eause higher pore loadings. 

Conclusions 

The adsorption of bovine heart eytoehrome c to 
the mesoporous siliea powder SBA-15 has been 
studied for different folding states of the protein 
and for a set of different ehemieal eonditions. 
The adsorption of folded eytoehrome c is usually 
explained in literature by a simple eleetrostatie in¬ 
teraetion model where the moleeules bind to the 
eharged siliea via eharged amino aeid residues 
on the protein’s surfaee. The paeking density is 
mainly defined by the repulsion between the in¬ 
dividual protein moleeules whieh is defined by 
their overall net eharge. By trying to extend this 
model to the adsorption behavior of unfolded ey¬ 
toehrome c, we find that this simple model fails to 
aeeount for the eomplieated dependenee of protein 
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binding on both the ionic strength and the pH of 
the buffer solution. We tried a different approach 
using the counterion release mechanism which 
readily explains the differences between the two 
folding states. However, the counterion release 
can not explain the residual adsorption observed 
for high ionic strength and the high pore loading 
at the protein’s isoelectric point. Only by com¬ 
bining electrostatic interaction, counterion release 
and van der Waals forces, we are able to explain 
the binding behavior of cytochrome c to silica. 
Apparently, protein binding to charged inorganic 
surfaces is caused by a complex interplay of differ¬ 
ent mechanisms. The application of these findings 
for the fractionation of mixed protein solutions 
will be part of a subsequent study. 
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Table 1: Fit parameters ni and a of the Langmuir-type isotherms in figure]^ and their respective margins 
of error. Folded samples are marked as/, unfolded as u. 



ni [jUmol/g] 

^nL [/imol/g] 

a [l//imol] 

Act [F/imol] 

pH 4.4/ 

3.37 

±0.40 

0.023 

± 0.006 

pH 4.4 u 

7.84 

±0.56 

0.036 

± 0.007 

pH 6.4/ 

22.9 

±7.4 

0.25 

±0.10 

pH 6.4 u 

9.00 

±0.22 

0.11 

±0.01 

pH 9.7/ 

25.9 

±4.3 

0.40 

±0.19 

pH 9.7 u 

12.8 

±0.9 

0.073 

±0.014 
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